A COMPUTER PROOF OF TURAN'S INEQUALITY 
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Abstract. We show how Turan's inequality P n (x) 2 — P n -\{x)P n +\(x) > for Legendre 
polynomials and related inequalities can be proven by means of a computer procedure. 
The use of this procedure simplifies the daily work with inequalities. For instance, we 
have found the stronger inequality \x\P n (x) 2 — P n —x (x)P n +i (x) > 0, — 1 < x < 1, 
effortlessly with the aid of our method. 



1. Introduction 
Turan showed in a 1946 letter to Szego that 

(1.1) A n (x) := P n (x) 2 - P n ^(x)P n+1 (x) > 0, x G [-1, 1], n > 1, 

where P n (x) denotes the n-th Legendre polynomial. Szego |)] gave four non-trivial proofs. 
Several authors have proven analogous statements for other families of orthogonal poly- 
nomials, and there is now a substantial body of literature [H| devoted to these and related 
results. The aim of the present note is to describe a computer algebra proof of Turan's in- 
equality that requires as input only the three term recurrence of the Legendre polynomials 
and the first two polynomials. Our method |I] is applicable to many other inequalities, 
including the following refinement of Turan's result, which appears to be new. 

Theorem 1.1. Let P n {x) denote the n-th Legendre polynomial. Then 

(1.2) \x\P n {xf - P n ^{x)P n+1 {x) > 0, x G [-1, 1], n > 1, 
with equality holding if and only if either x = and n is even, or \x\ = 1. 

2. The Proving Method 

We exemplify our proving method on the classical Turan inequality (|1.1|) . The idea 
underlying the method is complete induction on n. That is, we establish the induction 
step 

(2.1) An(x) > =>> A n+1 (x) > 0, x G [-1,1], n > 1, 

and afterwards we verify that the original inequality holds for n = 1. For proving (j2.1j) 
automatically, we construct a so-called Tarski formula whose truth implies the validity of 
the induction step. Tarski formulas are quantified formulas built via logical connectives 
from polynomial equations and inequalities over the reals. Upon replacing P n _i(x), P n (x), 
P n +i(x), P n+ 2(x) in ()2.1|) by indeterminates Y~_i, Y , Yi, Y 2 , we obtain the formula 

$:= (VY^YoMMeB.'.YZ-Y-M^O => Y 1 2 -Y Y 2 >0). 

Three things have to be remarked about this formula, (i) It can be decided algorithmi- 
cally whether or not a given Tarski formula is true. The classical decision procedure of 
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Tarski ^1] as well as the more efficient method of Cylindrical Algebraic Decomposition 
(CAD) due to Collins pQ are available for this purpose, (ii) If $ holds, then (j2.1j) is also 
true, for if the implication Y$ — Y-{Y\ > => Y 2 — Y0Y2 > holds for all real numbers, 
then it holds in particular for any real number P n+i (x) (n and x arbitrary) in place of Yj. 
(iii) Of course, $ is false. 

In order to make the proof go through, additional knowledge about the Legendre poly- 
nomials has to be encoded into the hypothesis part of formula Remarkably enough, 
in case of Turan's inequality it is sufficient to throw in the inequality's domain of validity 
and the classic recurrence ^0] 

(n + 2)P n+2 {x) = (2n + 3)xP n+1 {x) - (n + l)P n (x), n > 0, 

of the Legendre polynomials. This requires additional indeterminates N (representing n) 
and X (representing x). The refined formula is 

V N, X, F_i,y , Fi, Y 2 6 R : (N > 1 A -1 < X < 1 A 

(N + 2)Y 2 = -(N + 1)Y + (3X + 2NX)YlA 

(N + l)Y x = -NY_x + (X + 2NX)Y Q ) 

(Y 2 - Y^Y, > Yl - Y Y 2 > 0). 

Using CAD, this formula can be easily verified by the computer, and by the remarks 
above we may regard this as a computer proof for the fact that Turan's inequality holds 
for n + 1 whenever it holds for n. 

To complete the proof, we have to consider the induction base n = 1. Since Pq(x) = 1, 
P\{x) = x, and P 2 (x) = (3x 2 — l)/2, we just have to verify the obvious formula 

V X G R : -1 < X < 1 => |(1 - X 2 ) > 0, 

which we can again leave to the computer, if we want. 

Strict positivity of A n (x) for — 1 < x < 1 can be shown analogously. 

3. Remarks and Further Applications 

We have to dispel any hopes that our method yields a decision procedure for inequalities 
involving orthogonal polynomials or other special functions. Needless to say, there are 
many special functions that do not fit into our recursive framework. Roughly speaking, our 
procedure requires functions of n (and possibly other real parameters) such that the n-th 
value depends polynomially on a finite number, independent of n, of previous values. For 
instance, the Bernoulli polynomials B n (x) cannot be handled, since their recurrence "goes 
all the way back" . Even if an inequality is in the input class, our method may be doomed 
to failure because the sufficient condition that we check might not be satisfied although 
the conjectured inequality is true. In some cases the user can remedy this by inputting 
extra equations or inequalities that the functions in question satisfy. A third reason for 
failure are excessive computing time and memory overflows; this is what happened when 
we tried to reprove Gasper's extension [3j of Turan's inequality to Jacobi polynomials. 
Using Mathematica's implementation of CAD, we ran out of memory (3 GB) after having 
spent forty hours of CPU time (1.5 GHz). This is in contrast to the computation time 
needed for Turan's original inequality, whose proof was completed in just a second. The 
reason for this discrepancy are the additional two parameters appearing in the Jacobi 
polynomials. 

In view of the doubly exponential complexity of CAD, it is surprising that our method 
is able to verify quite a few inequalities from the literature with a reasonable amount 



of time. For instance, it is a matter of seconds to verify Turan's inequality also for the 
following quantities in place of P n {x): 

• Hermite polynomials H n (x) (for x £ R), 

• Laguerre polynomials L%(x) (for x > 0, a > 0), 

• normalized Laguerre polynomials L"(x)/L"(0) (for x > 0,a > —1), 

• differentiated Legendre polynomials -P^(x) (For —1 < x < 1; the inequality actu- 
ally holds for all igR, but our method fails outside [—1, 1]). 

None of these results are new. We can also prove the inequality 

71 — 1 

A„(z) > ^yA n _!(x), iG[-l ) l],n>2 ) 

which is due to Constantinescu 

Our method lends itself to playing around with conjectured inequalities; this is how 
Theorem 11.11 was obtained. Note that the absolute value function can be easily accom- 
modated by Tarski formulas. The cases where we claim equality in (jl.2j) follow from the 
well-known facts P n (l) = 1, P n {-1) = (-l) n , and 

JO, n odd, 

-Pn(O) = < (_i)n/2 , . 

These can also be proven automatically by the method described above. However, there 
are a lot of established methods available for which proving identities like these is offend- 
ingly trivial 0015! 

We believe that our method could become a helpful tool for researchers working with 
inequalities. It might not be capable of proving difficult inequalities that are of interest in 
their own right (Turan's inequality seems to be exceptional in this respect), but it might 
be helpful for proving elementary inequalities that appear as subproblems in the proof of 
more involved theorems. 
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